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Излагаются результаты вычислений частотного спектра радиосигнала, отраженного от турбулентного плазменного образования, 
создаваемого реактивным двигателем космического аппарата. Проводится сопоставление результатов численного моделирова- 
ния с известными экспериментальными данными по радиозондированию вслед факелу ракеты. 


Введение 

Изучение воздействия возмущений ионосфер- 
ной плазмы, создаваемых космическим аппаратом 
(КА), на распространение радиоволн началось с 
момента запуска первого искусственного спутника 
Земли как у нас в стране, так и за рубежом. Иссле- 
дованиям по проблеме воздействия на радиосигнал 
ракетно-ионосферных возмущений, образующих- 
ся в ближней зоне КА, посвящены работы [1-15]. В 
работе [1] приведен частотный спектр радиосигна- 
ла, отраженного от турбулентной струи выхлопных 
газов ракеты, движущейся в ионосфере, получен- 
ный для случая радиозондирования вслед ракете. 
Проведенное в указанной работе теоретическое ре- 
шение задачи о рассеянии радиосигнала на разре- 
женном облаке выхлопных газов в приближении 
однократного рассеяния не позволило объяснить 
полученные экспериментальные результаты. 

Следует отметить, что большое количество фак- 
торов, влияющих на структуру принимаемого ра- 
диосигнала (геофизические условия, рабочие час- 
тоты, ориентация радиотрасс, траектория движе- 
ния КА, вид топлива реактивного двигателя и т.д.), 
существенно затрудняет задачу однозначной инте- 
рпретации полученных экспериментальных дан- 
ных и не позволяет составить полную картину ис- 
следуемого явления. 

Цель данной работы состоит в исследовании 
воздействия ионосферных возмущений, создавае- 
мых реактивным двигателем в ближней зоне КА, на 
рассеяние радиоволн. 

1. Постановка задачи 

В работах [6, 7, 11, 15] численно исследовано 
прохождение электромагнитной волны через поток 
турбулентной слабоионизованной плазмы с неод- 
нородным профилем скорости и концентрации 
турбулентностей. На основе сопоставления резуль- 
татов расчетов с экспериментальными данными 
[16] в работе [15] сделан вывод о том, что турбуле- 
нтности сосредоточены в окрестности поверхности 
потока. В связи с этим представляет интерес поста- 
новка следующей задачи [4-7, 15]. 

Рассматривается рассеяние плоской радиовол- 
ны на турбулентных плазменных неоднородностях, 
сосредоточенных вблизи поверхности тел враще- 


ния: конуса, параболоида вращения и поверхнос- 
ти, образованной вращением кривой четвертого 
порядка [17]. 

Предполагается, что размеры тел вращения 
намного превосходят длину радиоволны и харак- 
терные размеры турбулентных неоднородностей, а 
час тота р адиоволны удовлетворяет соотношению 
гд е/ р - плазменная частота ионосферы 
на высоте движения КА, А, А 0І - возмущенная и 
невозмущенные концентрации заряженных частиц 
ионосферной плазмы. При этом коэффициент от- 
ражения от турбулентных плазменных неоднород- 
ностей близок к 1 [17]. 



Рис. 1 . Г еометрия задачи при радиозондировании внутрен- 
ней поверхности турбулентного тела вращения 


Полагается, что плоская электромагнитная вол- 
на падает в направлении, противоположном оси і 
тела вращения (конуса, параболоида или поверх- 
ности, образованной вращением кривой четвертого 
порядка), на поверхности которого хаотично распо- 
ложены турбулентные неоднородности. Они обла- 
дают как направленной скоростью Ѵ„, ориентиро- 
ванной вдоль образующей тела вращения, так и слу- 
чайной изотропной скоростью 8Ѵ, распределенной 
по нормальному закону. Решение задачи проводится 
в системе координат, связанной с КА, ее центр рас- 
положен в вершине тела вращения, а ось і направле- 
на вдоль его оси (рис. 1). В этом случае поле рассе- 
янной волны формируется в результате многократ- 
ного переотражения волны от внутренней поверх- 
ности турбулентного плазменного тела вращения. 

2. Стохастическое моделирование многократного 
рассеяния радиоволн на внутренней поверхности 
турбулентного плазменного тела вращения 

Существующие аналитические методы теории 
рассеяния волн в случайных дискретных средах 
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имеют ограниченную область применения [18, 19]. 
Поэтому для решения проблемы многократного 
рассеяния радиоволн на перемещающихся турбуле- 
нтных неоднородностях ниже применяется имита- 
ция вероятностной модели исследуемого явления. 
Турбулентным неоднородностям ставятся в соотве- 
тствие случайно независимые элементарные рассе- 
иватели. Полагается, что длина электромагнитной 
волны значительно меньше характерных размеров 
турбулентного плазменного образования. Предпо- 
лагаются выполненными условия применимости 
теории переноса излучения [18]. Для решения зада- 
чи используется метод Монте-Карло [20, 21]. 

Фронт плоской электромагнитной волны, пада- 
ющей на поверхность турбулентного плазменного 
образования, представляется источником синфаз- 
ных когерентных лучей. Наличие перемещающих- 
ся хаотически расположенных на поверхности тур- 
булентных неоднородностей приводит к тому, что 
фазы лучей после отражений от поверхности будут 
случайны. Энергия рассеянного сигнала в этом 
случае пропорциональна числу лучей, рассеянных 
в заданном элементе телесного угла. 

В плоскости фронта волны, падающей на тело 
вращения, выбирается луч, направление которого 
совпадает с направлением волнового вектора. На- 
чальные координаты луча распределены равномерно 
по раскрыву тела вращения в плоскости 1 = 7 ^, ограни- 
чивающей тело вращения. Луч продолжается до пе- 
ресечения с поверхностью тела вращения. Координа- 
ты первого пересечения луча с поверхностью опреде- 
ляются из условия распространения луча в направле- 
нии противоположном оси і, и выполнения вдоль его 
траектории соотношений х=сошІ и у=соші. 

Координаты второго и последующих пересече- 
ний луча с поверхностью находятся из совместного 
решения системы уравнений 


г=г 0 + \ Г-Щ, 

( 1 ) 

0(г)= 0 , 

( 2 ) 


где уравнение (2) является уравнением поверхнос- 
ти, на которой происходит рассеяние. В уравнении 
(1) известными являются % - координаты начала 
луча на рассеивающей поверхности и единичный 
вектор ё* в направлении его распространения. 

Ниже рассматривается рассеяние радиоволн на 
трех типах поверхности пл азмен ного образования: 
конической (^=р/1§0 о , р='$ 3 +/, Ѳ 0 - угол полура- 
створа конуса), в виде параболоида вращения 
(г=р 2 ) и поверхности, образованной вращением 
кривой, описываемой уравнением четвертого по- 
рядка (і=р 4 ). Величины ія р нормированы на зна- 
чение радиуса среза сопла г и и измеряются в отно- 
сительных единицах. 

При пересечении луча с поверхностью происхо- 
дит изменение направления его распространения в 
соответствии с задаваемой индикатрисой переизлу- 
чения турбулентных неоднородностей. Представ- 
ляет интерес рассмотреть диаграмму переизлучения 


рассеивателей квазизеркального типа, которая поз- 
воляет охватить широкий диапазон соотношений 
между характерным размером рассеивателей сі и 
длиной волны Я от ^/>>Я при узкой диаграмме пере- 
излучения рассеивателей до й«Х при широкой ди- 
аграмме, когда она стремится к изотропной. 

В этом случае единичный вектор в направлении 
распространения рассеянной волны определяется 
выражением 

е, = р - 2е ч (Ще д ), 
где 

% = {п + р)І\(п + р)\. 

Здесь п - нормаль к рассеивающей поверхности, 
- случайный вектор, касательный к поверхности в 
точке отражения. Выражение, используемое для мо- 
делирования вектора р*, записывается в виде 

р = (е„ С08 Ц/ + 8ІП у/ \ѵ) р фс~ { , 
где 

\Ѵ=(пхё„)/\пхё и [і; 

^-единичный вектор, направленный вдоль обра- 
зующей тела вращения; і//=2ях 2 - азимутальный 
угол в плоскости, касательной к поверхности в точ- 
ке рассеяния; х ь х 2 - случайные числа, равномерно 
распределённые в интервале от 0 до 1 ; /3 — коэффи- 
циент, характеризующий ширину диаграммы пере- 
излучения рассеивающей поверхности. В случае 
Р>> 1 элемент среды рассеивает изотропно. 

В процессе моделирования случайных траекто- 
рий лучей проводится вычисление и накопление в 
памяти компьютера величины доплеровского сдвига 
частоты, приобретенного волной в результате мно- 
гократного взаимодействия с перемещающимися 
турбулентными неоднородностями. При выходе лу- 
ча за пределы тела вращения запоминаются накоп- 
ленный им сдвиг частоты, направление его распро- 
странения и число отражений от поверхности. После 
многократного повторения процесса слежения за лу- 
чами и нормировки результатов получаются спект- 
ральные характеристики рассеянного сигнала [15]. 

На рис. 2 и 3 представлены частотные спектры 
сигнала, некогерентно рассеянного на внутренней 
поверхности тел вращения в случае квазизеркаль- 
ной и изотропной по полусфере диаграмм переиз- 
лучения турбулентных неоднородностей. По гори- 
зонтальной оси на рисунках отсчитывается величи- 
на безразмерного доплеровского сдвига частоты 
/•=(А/с)/(/ 0 Ѵ и ), где А/- сдвиг частоты, с - скорость 
распространения волны, - частота падающей 
волны, Ѵ и - скорость перемещения турбулентных 
неоднородностей вдоль образующей тела враще- 
ния. По вертикальной оси каждого спектрального 
кадра отсчитывается энергия сигнала, нормиро- 
ванная на максимум энергии в спектре. 

Значения параметров модели следующие: отно- 
шение максимального продольного размера плаз- 
менного образования к максимальному попереч- 
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ному равно % т /р т = 2; коэффициент зеркальности, 
характеризующий вероятность зеркального отра- 
жения лучей от поверхности а = 0,03; коэффици- 
ент Д характеризующий ширину квазизеркальной 
диаграммы переизлучения рассеивателей, прини- 
мал значение /3 = 1,5; относительная величина дис- 
персии скорости рассеивателей равна 5Ѵ/Ѵ= 0,2. 



Рис. 2. Частотный спектр радиосигнала, рассеянного на внут- 
ренней поверхности турбулентного тела вращения в 
случае квазизеркальной диаграммы переизлучения 

Первый столбец частотных спектров соответ- 
ствует рассеянию на конической поверхности, вто- 
рой - на параболоиде вращения и третий - на пове- 
рхности, образованной вращением кривой четвер- 
того порядка. В столбцах частотные спектры отли- 
чаются друг от друга значением полярного угла Ѳ от- 
носительно оси і. Накопление информации о доп- 
леровском сдвиге частоты проводилось в 10 ячейках 
по углу Ѳ, который изменяется от 0 до я/2. Нижние 
рисунки в столбцах соответствуют интегральному 
частотному спектру, полученному в результате сум- 
мирования всех лучей, вышедших из плазменного 
образования в диапазоне углов 0<Ѳ<л/2. 

Характерной особенностью спектра является 
наличие двух компонент: дискретной (при /=0), 
полученной в результате зеркального отражения 
радиоволны от поверхности, и сплошной, полу- 
ченной в результате рассеяния на перемещающих- 
ся рассеивателях. Справа относительно максимума 
сплошной составляющей наблюдается резкий спад 


до 0 энергии сигнала. Слева, с удалением от макси- 
мума, происходит плавное понижение энергии, 
причем ее значение всегда отлично от 0. Указанные 
особенности согласуются с результатами проведен- 
ного в [15] анализа величины доплеровского сдви- 
га частоты многократно рассеянного сигнала. 

Следует отметить, что с ростом порядка уравне- 
ния, описывающего поверхность, увеличивается 
значение энергии для отрицательного сдвига час- 
тоты /• (рис. 2, 3). Этот факт объясняется возраста- 
нием среднего числа актов рассеяния лучей на тур- 
булентностях, перемещающихся вдоль образую- 
щей плазменного образования, что приводит к уве- 
личению числа лучей с отрицательным значением 
величины доплеровского сдвига частоты. 
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Рис. 3. Частотный спектр радиосигнала, рассеянного на 
внутренней поверхности турбулентного тела враще- 
ния в случае изотропной диаграммы переизлучения 


Анализ результатов, представленных на рис. 2 и 
3, показывает смещение максимума сплошной сос- 
тавляющей спектра в сторону больших значений 
величины доплеровского сдвига частоты с ростом 
ширины квазизеркальной диаграммы переизлуче- 
ния рассеивателей. Дело в том, что с уменьшением 
коэффициента Д характеризующего ширину ква- 
зизеркальной диаграммы переизлучения рассеива- 
телей, средний вектор рассеяния радиоволны <? 
приближается к положению нормали к поверхнос- 
ти п, что и обеспечивает меньшую величину допле- 
ровского сдвига частоты на рис. 2 по сравнению с 
рис. 3 (поскольку (й’е^)=0) [15]. 
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3. Сопоставление результатов численного моделиро- 
вания с экспериментальными данными по радиозон- 
дированию вслед факелу ракеты 

В работе [1] приводятся результаты радиолока- 
ционных измерений в коротковолновом диапазоне 
частот при запуске ракеты. Радиозондирование 
проводилось вслед ракете, движущейся в верхних 
слоях атмосферы. Геометрия эксперимента приве- 
дена на рис. 4. Результаты измерений частотного 
спектра принимаемого радиосигнала представле- 
ны на рис. 5 сплошной кривой. По горизонтальной 
оси на рис. 5 отсчитывается величина безразмерно- 
го доплеровского сдвига частоты По оси ординат 
отсчитывается энергия сигнала, нормированная на 
максимум энергии в спектре. 

Контактная поверхность 



Рис. 4. Взаимное расположение ракеты и радиолокатора 
при измерениях [1] 

Как правило, частотный спектр радиосигнала, 
отраженного от ракеты с работающим реактивным 
двигателем, является бимодальным [1]. Появление 
одной из спектральных компонент обусловлено от- 
ражением радиоволны от тела самой ракеты, а ве- 
личина доплеровского сдвига частоты соответству- 
ет скорости движения ракеты [ 1 ] . 

Вторая спектральная компонента смещена от- 
носительно первой спектральной компоненты на 
величину доплеровского сдвига частоты и перек- 
рывает значительный участок спектра. Ее появле- 
ние обусловлено отражением радиоволны от пото- 
ка выхлопных газов реактивного двигателя. В этом 
случае радиоволна отражается от рассеивающих 
частиц, движущихся со скоростями, отличающи- 
мися от скорости ракеты-носителя примерно на 
величину предельной скорости истечения выхлоп- 
ных газов (~3 км/с) [1]. 

Проведенное в указанной работе решение задачи 
о рассеянии радиосигнала на разреженном облаке 
выхлопных газов в приближении однократного рас- 
сеяния не позволило объяснить полученные экспе- 
риментальные результаты. Дело в том, что в модели 
рассеяния радиоволн на турбулентном потоке, ана- 
лизируемой в работе [1], ширину частотного спект- 
ра рассеянного сигнала определяет разброс скорос- 
тей рассеивающих неоднородностей, который зна- 


чительно меньше предельной скорости истечения 
выхлопных газов. Поэтому указанной величины 
ширины частотного спектра рассеянного сигнала 
оказалось недостаточно для получения наблюдав- 
шегося в эксперименте (рис. 5) перекрытия диск- 
ретной и сплошной спектральной составляющих. 



Рис. 5. Сопоставление частотных спектров, полученных в ре- 
зультате расчетов, с экспериментально наблюдав- 
шимся частотным спектром, полученным при радио- 
зондировании вслед факелу ракеты [1] 

Результаты численного моделирования рассея- 
ния радиоволн на внутренней поверхности турбу- 
лентного плазменного тела вращения, проведенно- 
го в настоящей работе, приведены на рис. 5. Рас- 
сматривается случай радиозондирования вдоль оси 
полого турбулентного плазменного образования. 
На рис. 5 штриховой и штрихпунктирной кривыми 
представлены частотные спектры многократно 
рассеянного радиосигнала. 

Штриховая и штрихпунктирная кривые на 
рис. 5 получены при рассеянии радиоволн на пара- 
болоиде вращения (коэффициент зеркальности 
а=0,035) и поверхности, образованной вращением 
кривой, описываемой уравнением четвертого по- 
рядка (а=0,005). При этом параметры модели име- 
ли следующие значения: ъ т /р,„= 2; /3=1,5; 5Ѵ/Ѵ= 0,2. 
Полярный угол выхода лучей из плазменного обра- 
зования заключён в пределах от 0° до 9°. 

Экспериментально наблюдавшийся частотный 
спектр характеризуется совокупностью перекрываю- 
щихся дискретной (при /=0) и сплошной спектраль- 
ных составляющих. Расчёты показывают, что величи- 
на коэффициента /3, характеризующего ширину ди- 
аграммы переизлучения рассеивателей, существенно 
влияет на положение максимума сплошной состав- 
ляющей в спектре. При /3=1,5 положение указанного 
максимума на расчётных штриховой и штрих-пунк- 
тирных кривых близко к положению максимума на 
экспериментальной сплошной кривой. 

Коэффициент зеркальности а влияет на распре- 
деление энергии между дискретной и сплошной 
составляющими в спектре. Различным типам пове- 
рхностей в данном случае соответствуют такие зна- 
чения коэффициента а, при которых соотношение 
между величинами дискретной и максимума 
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сплошной составляющей близко к соответствую- 
щему соотношению на экспериментальной кривой. 

Подобие формы экспериментального спектра и 
кривых, полученных в результате численного моде- 
лирования, позволяет сделать вывод о том, что та- 
кой спектр у отражённого сигнала появляется в ре- 
зультате многократного рассеяния электромагнит- 
ной волны на внутренней поверхности турбулент- 
ного плазменного тела вращения. При этом запол- 
нение провала между максимумами наблюдавше- 
гося в эксперименте бимодального частотного 
спектра объясняется накоплением величины доп- 
леровского сдвига частоты радиоволны в процессе 
многократного взаимодействия с турбулентными 
неоднородностями, перемещающимися вдоль об- 
разующей плазменного тела вращения. 

Заключение 

1. Методом Монте-Карло исследовано многок- 
ратное рассеяние радиоволн на внутренней по- 
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